Baryon Resonances in a Chiral Hadronic Model for the QCD Equation of State 
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In this paper we study the influence of hadronic resonances on the phase diagram calculated 
with an effective chiral flavour SU(3) model. We show that varying the couplings of the baryonic 
resonances to the attractive scalar and the repulsive vector flelds has a major impact on the order 
and location of the chiral phase transition and the possible existence of a critical end point as well 
as on the thermodynamic properties of the model. Furthermore, we study (strange) quark number 
fluctuations and show the related susceptibilities both at zero baryochemical potential and when 
crossing the phase transition line at three different points in the T-^ plane. We obtain the best 
agreement with current lattice data if we choose a rather strong vector coupling which in our model 
limits the phase transition to a smooth crossover and implies the non-existence of a critical end 
point. 



I. INTRODUCTION 

Currently there are many efforts both on the theoreti- 
cal, as well as on the experimental side to gain knowledge 
about the phase diagram of strongly interacting nuclear 
matter. Experiments with highly energetic colliding gold 
nuclei at the Relativistic Heavy Ion Collider (RHIC) sug- 
gest that a new state of matter, comparable to a nearly 
perfect fluid, is created at high temperatures [H-ll]. While 
these experiments at RHIC and the Large Hadron Col- 
lider (LHC), performed at very high beam energies, aim 
at the high-temperature, low-density region of the phase 
diagram, there are also experiments in which higher bary- 
onic densities may be reached, as for example planned at 
the upcoming Facility for Antiproton an Ion Research 
(FAIR) at GSI. Common to all these experiments is the 
search for the phase transition of strongly interacting 
matter from a confined state of hadrons at low tempera- 
tures to a state with deconfined quarks and gluons (QGP) 
as it is predicted for high densities and high tempera- 
tures [H. This phase transition is commonly referred to 
as deconfinement transition [6| • The symmetries of quan- 
tum chromodynamics (QCD), however, imply another, 
chiral phase transition at which the chiral symmetry is 
restored and the masses of the baryons, or constituent 
quarks, vanish 043- For this transition, on which we 
will mainly focus our studies, the chiral condensate a 
acts as the corresponding order parameter. 

Since QCD can not be treated perturbatively at low 
temperatures, most information we have about the phase 
diagram of QCD matter today comes from lattice QCD 
calculations |lC)| - [l9| . In the region of vanishing bary- 
ochemical potentials ^ = and finite temperature, lattice 
QCD yields reliable results that show a smooth crossover 
transition [l^. For finite potentials /i > 0, however, re- 
sults from lattice QCD are not available because of the 
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failing of standard Monte Carlo sampling methods due 
to the so-called sign-problem. Currently there are vari- 
ous methods to extend the lattice results from a = to 
the region of small chemical potentials [ll|, [12, [20l - [29| . 
Some lattice QCD groups suggest that the phase transi- 
tion becomes first order at a critical endpoint J^, .301- 
[32j with its coordinates varying considerably for different 
calculations. Other more recent results [H] favour the ex- 
clusive existence of a crossover transition for a wide range 
of baryon densities, implying the non-existence of such a 
critical point. Moreover, the exact position of the phase 
transition in the phase diagram is subject of a lively and 
ongoing debate. 

Other common theoretical approaches to QCD are ef- 
fective models to study specific properties of QCD mat- 
ter 0, |33 - |40| . In our approach we will use an effective 
model for the QCD equation of state to study the prop- 
erties of the phase diagram of nuclear matter. Our model 
is able to reproduce the well known saturation properties 
of nuclear matter [4l|, . 

The degrees of freedom in our model are all known 
baryons and baryonic resonances up to masses of m = 
2.6 GeV. This approach resembles an interacting hadron 
resonance gas (HRG). The HRG is known to give a good 
description of thermodynamic quantities at low temper- 
atures (T < Tc) and was often used to reproduce lat- 
tice QCD results for thermodynamics and fluctuations of 
conserved charges at low temperatures [H, lisl - lisj . The 
particle production in heavy ion collisions may also be 
reasonably well described with use of the HRG [49 'ss']. 
Therefore, our first approach to use only hadronic degrees 
of freedom seems to be appropriate up to at least tem- 
peratures in the range of T^. For higher temperatures, at 
the latest when the deconfinement phase transition sets 
in, this approach breaks down and, in a future step, par- 
tonic degrees of freedom need to be taken into account 
as it was done in Refs. [13, [Ell for example. 

This paper is organised as follows. First we introduce 
the chiral SU{Z) model in section II. In section HI, we 
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present results for the order parameter of the chiral con- 
densate a at zero and nonzero baryochemical potentials, 
thermodynamic quantities from the model, and quark 
susceptibilities at different points in the phase diagram. 
Susceptibilities are interesting quantities since they re- 
semble fluctuations of conserved charges which itself are 
closely linked to phase transitions. They offer an effec- 
tive possibility of comparison between results from the- 
ory and experiment because the susceptibilities can be 
related to measured fluctuations of particle production. 
Of particular interest are susceptibility ratios since they 
are not dependent on the volume or the impact parame- 
ter of the underlying system [56| . This work closes with 
a conclusion in section IV. 



II. MODEL 

In our model, a S'J7(3)-flavour sigma-omega model us- 
ing the non-linear realization of chiral symmetry (see 
Refs. [13, EH, m, [131 for a comprehensive review), the 
Lagrangian in mean field approximation has the form 



(1) 



Here, the first term represents the kinetic energy of the 
hadrons, the terms 



(2) 



describe the interaction of the baryons with the scalar 
mesons tr, ( (attractive interaction, see Eq. ([5])) and the 
vector mesons w, (j) (repulsive interaction) respectively. 
The summation index i runs over the baryon octet {N, A, 
E, S), the baryon decuplet (A, E*, S*, f2), and aU heavier 
resonance states up to masses of tojv* = 2600 MeV. We 
only include hadronic resonances whose existence is con- 
sidered to be v ery likely according to the Particle Data 
Group listings [58| where they are recorded with a min- 
imum three-star rating. Since the listed masses of some 
heavy resonances may cover a broad range, all particles 
are included with their average mass. 
The third term of the Lagrangian 
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(3) 



V 4 ^/2 %/2 



(4) 



fc2 + C + fcsfr'C + ki In 
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includes the self interactions of the vector mesons and the 
scalar mesons together with the last two terms describing 
the explicit symmetry breaking. 
The effective masses of the baryons 



(5) 



are created by the coupling of the baryons to the scalar 
meson fields (i.e. the non-strange chiral condensate a and 
its strange equivalent (), together with an explicit mass 
of at least SniN = 150 MeV. In this way at high temper- 
atures and baryonic densities the decreasing u-field leads 
to smaller baryon masses and thus to the restoration of 
chiral symmetry. Thereby the nucleons which have the 
smallest explicit mass have lost roughly 45% of their vac- 
uum mass at Tc- Since our model only includes hadronic 
degrees of freedom it can only be applied in the hadronic 
regime up to temperatures slightly above Tc- This also 
prevents baryonic masses from getting to small so that 
the mean field approximation should remain valid and no 
fiuctuations need necessarily to be taken into account as 
it was for example done in Refs. [sol. l60j. 

The effective masses of the pseudoscalar mesons and 
vector mesons are given by the second derivative of the 
mesonic potential Vmcson = —^vcc ~ ~ ^esb with re- 
spect to the respective mesons at the minimum of the 
grand canonical potential (see Eq. ©) [6l| 



(6) 



£.3 = '^,V,v\K,K; p,uj,if,K*,K*. 



The couplings of the baryon octet to the mesonic fields 
and the mesonic potential are chosen in such a way as to 
reproduce the well-known vacuum masses, the nuclear 
ground state properties (e.g. the correct binding energy), 
and the asymmetry energy. The coupling strengths of 
the baryons of the decuplet and all heavier resonances 
are scaled to the nucleon couplings via the parameters 
Ts, according to 



(7) 
(8) 



In this paper we systematically analyse the influence of 
the baryonic resonances on the hadronic matter proper- 
ties. Therefor we adjust the vector coupling parameter 
r„ , controlling the abundance of the baryonic resonances 
at finite baryochemical potentials in the model, in order 
to study the infiuence on the resulting phase diagram 
and the thermodynamic properties of the model. We 
use this one-parameter approximation in order not to be 
swamped by a plethora of of unknown coupling constants 
of the various hadronic multiplets. In this work we set 
the scalar coupling parameter fixed to = 0.97 which 
ensures a smooth crossover phase transition at zero bary- 
ochemical potential (see Fig. [T] (b) for a study of the ef- 
fect of Ts). In general the scalar couplings are fixed by 
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reproducing the particles' vacuum masses (except for the 
explicit mass term Srrii). 

As a reference we also perform calculations for an ideal 
HRG that is not interacting with the mesonic fields. In 
this particular case all baryon couplings to the fields gB 
are set to zero and the masses of all particles arc fixed at 
their tabulated vacuum expectation value. 

The grand canonical potential of our model takes the 
form 



th 



(9) 



with the thermal contribution of the hadrons in the model 



ieB 



+ ln[l + e-T(BrW+M*)]^ (10) 



-T 



^ (27r)3 



J <fk In [l-e-^(^^*('=)-^^ )" 



Here the sums run over all baryons B (anti baryons are 
explicitly included in the second term of the first inte- 
gral) and all mesons M in the model, 7^ j- stands for the 
spin-isospin-degeneracy factor of the respective particle 

species i,j and E, 

tide energies 



j j(fc) = y fc2 + TO*^ for the single par- 
The effective baryochemical potential is 



defined as ^* 



Qi^cj). The grand canonical 



potential then leads to the thermodynamic quantities of 
the system, i.e. the pressure p and the energy and entropy 
density e, s, together with the densities of the particular 
particle species pi. 

Another effect we include in our model (following the 
works of Refs. [H, H^) is the effective suppression of 
baryonic states at high densities due to excluded volume 
effects (63-66] . While so far all particles were regarded 
as point like, we now can introduce for every single par- 
ticle j an average finite volume v^^ excluded from the 
total volume of the system V . In a first and very basic 
consideration this excluded volume in a non-relativistic 
definition takes the form 



V; 



^ -^l7r(2r)^ 
lai 6 



(11) 



where r is the mean radius of all particles in the model 
and the parameter defines the excluded volume for the 
different particle species i. In our case we set ab.m — 1 
as a first approximation for all baryons and mesons. This 
modification leads to the altered chemical potential 



(12) 



where P is the sum over all partial pressures. All ther- 
modynamic observables must now be expressed in terms 
of T and /i^-. Furthermore, a volume correction factor for 
each particle species 



V 
V 



Pi 



(13) 



with VI the volume not being occupied, has to be intro- 
duced in order to express the densities in a thermody- 
namically consistent way 



P't = ftPt, 

i 



(14) 
(15) 

(16) 



The inclusion of the excluded volume effects also has an 
impact on the nuclear ground state properties and thus 
this approach is considered only as a first attempt to 
study the impact of the suppression of baryons on the 
phase transition and on the thermodynamic properties 
of the model. However, reasonable nuclear ground state 
properties can be achieved by recalibration of the param- 
eters as it was shown in Ref. 1671. 



III. RESULTS 

A. Order Parameters 

We start our investigation of the model properties with 
the calculation of the order parameter for the chiral 
transition a at zero baryochemical potential pB- Fig- 
ure [T] (a) shows the normalised chiral order parame- 
ter as a function of the temperature together with data 
from lattice QCD calculations. Here and in the follow- 
ing, the lattice data we refer to were obtained by differ- 
ent collaborations using various lattice actions (asqtad, 
hisq, p4 and stout) and t emp oral spacings of the lattice 

(iv^ = 4 - 12) p, m, m, nHzi- 

With the scalar couplings fixed to = 0.97 we ob- 
tain a smooth crossover in a for both the model with 
excluded volume (green dashed line) effects due to fi- 
nite size of the particles and without it (red solid line). 
The critical temperature, defined as the point with the 
largest increase in a, is found to be Tc = 164 MeV with- 
out the excluded volume effect and Tc = 174 MeV if 
the finite size effects are taken into account. The results 
for a without the excluded volume effects are in quali- 
tatively good agreement with lattice QCD calculations 
which predict a smooth crossover at /is =0 and crit- 
ical temperatures in the range from Tc = 155 MeV to 
200 MeV [H, m, m, Hi-lzl [73, where the newest con- 
tinuum extrapolated data from the Hot-QCD and the 
Wuppertal-Budapest Collaborations predict critical tem- 
peratures close to r = 160 MeV consistently. For tem- 
peratures below Tc the slope of the chiral condensate 
from our model is relatively steep and deviates from those 
suggested by lattice QCD calculations. This is mainly 
due to neglecting contributions from pseudoscalar me- 
son self interactions which are important at low temper- 
atures [ill. The slope of a could be leveled by including 
the TT self interaction in the model [63, u3 . 
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Figure 1. (Color online) (a) Normalised order parameter for 
the chiral condensate a/ao as a function of T at /is = 
with (green dotted line) and without (red solid line) excluded 
volume corrections compared to lattice data. The blue dotted 
line at a/ao = 1 shows the reference line of the non- interacting 
ideal hadron resonance gas. Here and in the following plots 
the lattice data are shown as calculated with the asqtad [6g| . 
the hisq [H.ItO], the p4 [H, [H, [tJI , and the stout action^, 
[zl: [z3 on lattices with different temporal extent A'^^ . Panel 
(b) shows a/ao for different strength of the scalar coupling 
together with results from our model without any resonance 
states. 



The impact of the heavy resonances states can be seen 
from Fig. [T] (b) where the black dots show the normalised 
chiral order parameter as a function of the temperature as 
calculated with our original model which does not include 
any hadronic resonance states [4l|, It shows that 

without the hadronic resonance states the transition is 
much smoother and exhibits a higher critical temperature 
close to Tc = 200 MeV. The same effect is achieved by 
decreasing the strength of the attractive coupling of the 
baryons to the scalar cr-field. This behavior is also shown 
in Fig. [T] (b) for values of rg from to 0.9. The weaker 
the scalar coupling is, the flatter is the transition and 
the critical temperature moves to higher temperatures. 
For vanishing scalar couplings (red line) the results are 
the same as in the original model. Note that for small 
values of Vs one would have to introduce large values 
for the explicit mass term as the sum of the meson-field 
generated mass and the explicit term have to reproduce 
the vacuum masses of the states. 

Next we extend this study of a to nonzero baryochem- 
ical potentials, i.e. the whole T-^b plane. One of the 
major interests in the study of the phase diagram of 
strongly interacting matter is certainly to obtain infor- 
mation about the chiral and deconfinement phase transi- 
tions and the search for a possible critical end point which 
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Figure 2. Phase transition lines of the chiral order parameter 
a in the T-fiB plane for different values of r„ together with 
the nuclear ground state at /is = 939 MeV, T = and the 
liquid gas first order phase transition. The point on each 
phase transition line depicts the critical end point, separating 
a first order phase transition (solid lines) from a crossover 
(doted lines). Note that for values Vv > 0.6 there is solely a 
very broad crossover in the whole T-^b plane. 



divides the crossover phase transition at vanishing bary- 
ochemical potential from a first order phase transition at 
finite chemical potentials. Since lattice QCD calculations 
are systematically limited to /is = due to the so-called 
sign problem, different lattice QCD groups use various 
methods to extend their results to non-vanishing poten- 
tials. Unfortunately, this did not lead to a consistent pic- 
ture for the critical end point until now. While results 
from (2+l)-flavour QCD calculations suggested the criti- 
cal end point at a critical quark chemical potential in the 
range from /ig" w 725 MeV :2(1] to /i^"* w 140 MeV [zl 

(with 3/iq = Ms) ^-iid Mb'* ~ 360 MeV more recent 
studies of the Wuppertal-Budapest Collaboration [s^ 
doubt the existence of a critical end point at all and sug- 
gest a broad crossover phase transition over the whole 
T-pB plane. 

In Fig. [5] we show the results from our model with- 
out excluded volume effects for different values of the 
coupling constant of the repulsive vector interaction ry. 
In this figure the lines represent the biggest increase of 
the (T-field when varying T and fi. While a discontin- 
uous first order phase transition in a is depicted by a 
solid line, a smooth crossover in a is drawn with a dot- 
ted line. The point on each phase transition line stands 
for the critical end point for the specific choice of r„. We 
find that increasing the vector coupling strength leads to 
a stronger suppression of heavier particles and thus to 
a phase transition at higher chemical potentials. Note 
that for all coupling strengths smaller than r„ = 0.4 the 
nuclear ground state a,t fiB — 939 MeV, T = and the 
first order liquid gas phase transition are located above 
the phase boundary in the chirally restored phase. There- 
fore, within our simple model approximations those small 
vector couplings do not lead to physically reasonable so- 
lutions. Thus, we conclude according to our model the 
smallest critical chemical potential (for r„ = 0.4) is lo- 
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Figure 3. (Color online) Energy density (a) and the pressure 
(b) divided by T** at = as a function of the temper- 
ature. Depicted are results for the interacting hadronic gas 
with (green dashed line) and without (red solid line) excluded 
volume corrections and for the non-interacting ideal hadron 
resonance gas (blue dotted line) together with lattice data. 



cated at /i'^'' « 840 MeV which is two times higher than 
the value suggested in Ref. [73l |. 

For vector couplings of = 0.6, in our model, there 
is a very small temperature range of first order phase 
transition going up to T"^*^ « 20 MeV, at higher temper- 
atures there is a broad-ranged crossover transition. For 
all higher vector couplings, the first order phase transi- 
tions vanish completely and only a broad crossover region 
remains. 



B. Thermodynamics 

The thermodynamic quantities, as predicted by our 
model at fis = are depicted in Fig. [31 Panel (a) shows 
the energy density and panel (b) shows the pressure both 
divided by T**, in comparison to lattice data. Here and 
in the following plots, the solid red curve represents the 
results from our regular model including all baryonic res- 
onances, the green dashed curve depicts the results from 
the model with the excluded finite size volume effects 
taken into account and the blue dashed line stands for 
the non-interacting ideal HRG without finite size effects. 

For the standard case in our model, at Tc the energy 
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Figure 4. (Color online) The interaction measure, three times 
the pressure subtracted from the energy density over T*, as 
a function of the temperature at fiB ~ for the interacting 
hadronic gas with (green dashed line) and without (red solid 
line) excluded volume corrections and for the non-interacting 
ideal hadron resonance gas (blue dotted line) compared to 
lattice data [H, [H, [H, M, M, [ill, Izl, S [zl- The data for 
the hisq action is shown for two different values of the light 
quark mass mi = 0.05 ms and 0.20 ms, where rus stands for 
the strange quarks mass. 



density and the pressure rise rapidly due to the emerg- 
ing abundance of particles in the system because of their 
decreasing masses. Besides a slightly higher critical tem- 
perature, the suppression of heavier particles caused by 
the finite size effects leads to a much smaller maximum 
of the energy density at Tc and a decrease for higher 
temperatures. In this case, three times the pressure di- 
vided by exhibits only a slight increase with a maxi- 
mum around Tc and decreases smoothly again for higher 
temperatures. The HRG results rise monotonically, as 
expected. 

Figure |4] shows the so-called interaction measure de- 
fined as the energy density minus three times the pres- 
sure divided by as a function of the temperature. 
Again, the results from our regular model (red curve) 
show a rapid increase of this quantity at Tc, the HRG 
rises monotonously being in qualitatively good agreement 
in the re gion around the critical temperature with lattice 
data [isl. IgsI, [Tlj which shows a peak slightly above Tc- 

Comparing our results for the thermodynamic quan- 
tities to lattice data, we state that a good agreement 
up to the critical temperature is achieved for our results 
from the interacting HRG without excluded volume ef- 
fects and the continuum extrapolated stout action with 
physical quark masses. However, for temperatures even 
below the critical temperature, we overshoot the results 
from all other lattice actions. 
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Figure 5. (Color online) Second order (a) and fourth order 
(b) quark number susceptibilities as a function of T/Tc at 
/xs = 0. Results for the interacting hadronic gas {rv = 0.8) 
with (green dashed line) and without (red solid line) excluded 
volume corrections and for the non-interacting ideal hadron 
resonance gas (blue dotted line) are shown compared to lattice 
data HEil. 



C. Susceptibilities 

As mentioned above, lattice QCD calculations are 
largely restricted to zero chemical potential because of 
the sign problem. However, thermodynamic quantities 
at non-vanishing potentials may be obtained by Taylor 
expansion methods as described in Refs. [2^ lisl [tsI. [79| . 

The pressure at a specific point in the phase dia- 
gram p(T,fiB) with preferably small chemical potentials 
is calculated by expanding the pressure at p{T, fj,B = 0) 
around ^lb/T 



P(T, hb) 

rp4 



n=0 



T I 



(17) 



with the Taylor coefficients 



c„(T) 



1 9"(p(T,a.b)/T4 



n\ d{iiBlTr 



(18) 



3=0 



The coefficients Cn{T) are proportional to the quark 
number susceptibilities reflecting quark number fluctu- 
ations [131 ■ 

Similarly, by expanding the pressure with respect to 
the strange chemical potential fig, the first order strange 
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Figure 6. (Color online) Ratio of the fourth order to the sec- 
ond order quark number susceptibility as a function of T/Tc 
at /is = for the three different cases as described in Fig. [5] 
together with lattice data [H, H, [tJ, . 



1.3 1.4 quark susceptibility Xs (T) can be found to be 



XsiT) = T' 



dpi 



(19) 



In contrast to lattice QCD, with our model we are able to 
calculate the pressure at any point in the phase diagram 
and can therefore numerically calculate the correspond- 
ing susceptibilities. 

In Fig. Owe show the second order C2 (a) and fourth 
order C4 (b) quark number coefficients in comparison to 
lattice QCD data [H, H, ItTI, III] at /xs = 0. Both, 
the model with and without the excluded volume effects 
show a maximum of C2 at the critical temperature and 
a smooth decrease for higher temperatures. The fourth 
order coefficient C4 exhibits a narrow peak at Tc- As ex- 
pected, the ideal non-interacting HRG susceptibilities do 
not show any peak due to the absence of a phase transi- 
tion. The susceptibilities of the ideal HRG come closest 
to the lattice data with the p4-action. However, a com- 
parison to the not yet available susceptibilities from the 
continuum extrapolated stout action would be very in- 
teresting since the lattice results with the stout action 
gives the best agreement for the thermodynamic quanti- 
ties and the strange quark susceptibility (Fig. [7]) . 

In our study, we found that the strength of the repul- 
sive vector coupling has a major influence on the extent 
of the fluctuations at the phase transition, even though 
the derivative is performed at fiB =0. This effect was al- 
ready pointed out in Refs. [Ull^. A lower vector repul- 
sion leads to significantly higher fluctuations of conserved 
charges at Tc- In our study we set the vector coupling to 

— 0.8 in order to match our results with lattice data 
up to T « Tc. This vector coupling strength stands in 
contradiction to the one we found earlier to be suitable 
to reproduce a critical end point in the region suggested 
by lattice QCD results (r„ < 0.2, cf. Fig.©. 

In Fig. ini we show the ratio of the fourth order to the 
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Figure 7. (Color online) Strange quark number susceptibility 
divided hy at fj,B = as a function of T for the three 
different cases as described in Fig. \5\ together with lattice 
data m, 0, m, [tO, [t^ ■ The vector coupling strength is set 
to r„ = 0. 



second order baryon number susceptibility. For the reg- 
ular model without finite size corrections the ratio is in 
line with the completely flat curve of the HRG up to the 
critical temperature. The curve exhibits a discontinuity 
at Tc with a sudden decrease above. For higher temper- 
atures the curve is mostly constant around C4/C2 ~ 0.3. 
Due to the effective suppression of degrees of freedom at 
the phase transition, the results with excluded volume 
effects have a significantly different shape with a much 
smoother curve showing only a minimum peak at Tc- 

The strange quark number susceptibility Xs divided 
by is shown in Fig. [71 Again, the three cases of 
our model yield significantly different results. The curve 
for the non-interacting ideal HRG rises again monoton- 
ically as expected. The regular model with all interac- 
tions switched on shows a massive rise at . The results 
from our interacting model is once again in good agree- 
ment with results from lattice QCD with the continuum 
extrapolated stout action up to temperatures of Tc. In 
the case with the excluded volume effects, in particular 
heavy strange particles are suppressed and thus the curve 
is much lower in this case showing a maximum at Tc and 
a slow decrease for higher temperatures. 

In a next step we study the susceptibility coefficients at 
non-zero baryochemical potentials. Therefore, we com- 
pare the coefficients in three different regions of the phase 
diagram, namely in the crossover region, at the critical 
end point, and at the first order phase transition. For this 
purpose, in a first attempt we set the vector coupling to 
Tv = 0, because this leads to a clearly observable and 
well defined phase transition line and a critical end point 
at /i^'* « 220 MeV; later on we will also calculate sus- 
ceptibilities with more realistic and stronger vector cou- 
plings. Note, that for higher values of r„ this subdivision 
of the phase diagram into three clearly distinguishable 
regions no longer applies. For a vector coupling strength 
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Figure 8. (Color online) Second order quark number suscep- 
tibilities C2 for vector couplings r„ = (a) and ri, = 0.8 (b). 
The quantities are calculated going perpendicularly through 
the phase transition line at three different points with /^s 7^ 0. 
For rv — the blue line depicts the crossing of the phase 
boundary in the first order phase transition regime, the green 
line in the crossover regime, and the red line shows C2 going 
directly through the critical end point (see text for more in- 
formation). Shown is also the reference of the ideal hadron 
resonan ce g as at the critical end point and lattice data at 



of Ty > 0.6, our model exhibits only a broad crossover 
transition in the whole T-fiB plane (cf. Fig. [2]). 

With our model we then calculate the second and 
fourth order susceptibilities along straight lines perpen- 
dicular to the phase boundary going through the points 

on the phase transition line ^^■^ = 29.3 MeV, T^^^ = 

(2) 

161.7 MeV for the crossover region, /i^' = 216.0 MeV, 
T(2) = 150 MeV at the critical end point, and /i^^ = 
489.0 MeV, T^^) = 105 MeV for the first order phase 
transition region. 

In Fig. [8] we show the second order coefficients for the 
three different regimes for = in panel (a) and for r„ = 
0.9 in panel (b), together with lattice data [S. lesl. [tiI. [73| 
calculated at vanishing baryochemical potential /is = 0. 
Note, that choosing a vanishing vector coupling of the 
baryon resonances leads to unreasonably high suscepti- 
bilities at Tc and must be regarded as a limiting test case. 
The green dashed curve shows the susceptibilities in the 
crossover region at low hb- Due to the smooth transition 
of the quantities, the susceptibility only shows a small 
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Figure 9. (Color online) Ratios of the fourth order to the sec- 
ond order quark number susceptibilities for vector couplings 
— (a) and r.y — 0.8 (b) at different regimes of the phase 
diagram. The line styles are as described in the previous fig- 
ure. Please note the logarithmic scale in the upper panel (a). 



maximum (rise) at Tc for — 0.9 (r^ — 0). The suscep- 
tibilities at the discontinuous first order phase transition 
(blue dashed line) show the highest maximum for both 
values of expected. At the critical end point (red 

solid line), the maximum of C2 is located in between the 
results from the crossover and the first order phase transi- 
tion region. As shown previously, the susceptibilities for 
the non-interacting ideal HRG (pink dashed line) show 
a monotonous rising behaviour not being affected by a 
varied vector coupling strength. 

The ratio of the fourth order to the second order coeffi- 
cients in the three different regions of the phase diagram 
is shown in Fig.[S]for vector couplings r^, = in panel (a) 
and — 0.9 in panel (b). Once again, the height of the 
peak at the critical temperature shows a major depen- 
dence on the strength of the vector coupling. Comparing 
the curves of the three different regimes of the phase di- 
agram, the susceptibilities show a significantly different 
slope around Tc. At non- vanishing baryochemical poten- 
tials the ratios deviate significantly from the HRG results 
what is similar to the findings of suppressed fluctuations 
in higher order cumulants in the vicinity of the transition 
region at zero baryochemical potential [s^l. In our case 
the fluctuations are suppressed by the repulsive vector 
interactions. This fact could help to probe the location 
of the phase transition and a possible critical end point 



experimentally by extraction of the susceptibilities from 
data. 

However, the good agreement of our results for the 
susceptibilities in the crossover region given a vector cou- 
pling strength of = 0.9 with the lattice data for vanish- 
ing baryochemical potential (green dashed line in Fig. [3] 
(b)) together with the results for the chiral order param- 
eter cr in the whole T~fi plane (Fig. [5]), leads us to the 
conclusion, that the existence of a critical end point ac- 
cording to our model is very questionable. The results 
from our model suggest that for reasonable values of the 
vector coupling, throughout the entire phase diagram the 
phase transition is a smooth crossover, ruling out the ex- 
istence of a critical end point. This finding corresponds 
well with the results from most recent lattice data 13311. 



IV. CONCLUSIONS 

We presented an effective chiral SU{3) model for the 
QCD equations of state. In this sigma-omega model we 
included all known hadrons up to resonances with masses 
of 2.6 GeV together with a parameter r„ that controls the 
coupling strength of the baryons to the repulsive vector 
meson field. Furthermore, we include a finite size effect 
that effectively suppresses heavy particles at higher den- 
sities. 

Using this model, at zero baryochemical potential 
we found a smooth crossover phase transition for the or- 
der parameter of the chiral condensate a with a criti- 
cal temperature in the range from Tc = 164 MeV to 
Tc = 174 MeV depending on the excluded volume ef- 
fects being taken into account. These results are in good 
agreement with various data from lattice QCD. Extend- 
ing this study to finite /i^, we show the strong depen- 
dence of the phase transition, i.e. position and order, on 
the vector coupling strength. For reasonable values of r„ , 
we find that the phase transition is a smooth crossover 
in the whole T-fi b plane and that there is no critical end 
point. 

We also show the thermodynamic quantities from the 
model and calculate the quark number and strange quark 
number susceptibility coefficients at different values of 
fi B ■ The susceptibilities show a good qualitatively agree- 
ment with lattice data if a sufficiently strong vector cou- 
pling is chosen. This finding underlines the model sug- 
gestion of the non-existence of the critical end point. If 
the susceptibilities are extracted at different positions on 
the phase boundary, we show that their significantly dif- 
ferent behavior may be used to distinguish the order of 
the phase transition at a given point. 
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